Electron loss of fast projectiles in the collisions with molecules 
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The single and multiple electron loss of fast highly charged projectiles in the collisions with neutral 
molecules are studied within the framework of a nonperturbative approach. The cross sections for 
single, double, and triple electron losses are calculated for the collision system Fe q+ — > N2 (g=24, 25, 
26) at the collision energies 10, 100, and 1000 MeV/u. The effects caused by the collision multiplicity 
and the orientation of the axis of target molecule are treated. It is shown that collision multiplicity 
effect leads to considerable differences for the cases of perpendicular and parallel orientations of the 
molecular axes with respect to the direction of the projectile motion, while for chaotic orientation 
such effect is negligible. 
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I. INTRODUCTION 

Stripping, or electron loss by heavy projectiles in the 
collisions with atoms has been subject of extensive stud- 
ies during last decade. Various approach for the theo- 
retical study of the projectile electron loss cross sections 
in the collisions with neutral targets are presented in the 
Ref. \M and in the monograph by A. Voitkiv and J. Ull- 
rich [2[ where detailed discussion of the experimental re- 
sults can be found. In theoretical study of such processes 
it is important the treatment of the effects caused by 
strong field of the projectile and accounting the transi- 
tions in the target state. Such a treatment requires us- 
ing nonperturbative methods for the calculation of strip- 
ping cross sections [3, 4]. Recently such an approach was 
developed on the basis of sudden perturbation approxi- 
mation and used for the study of electron losses by fast 
highly charged projectiles with neutral gas atoms [a, El- 

Despite the fact that considerable progress is made in 
theoretical study of electron loss by fast projectiles, most 
of the treatments are restricted by considering atomic 
targets. 

However, study of the electron loss by fast highly 
charged projectile in the collisions with molecular tar- 
gets is of fundamental and practical importance because 
of the number of new effects that cannot be observed in 
the case of atomic targets. 

One of such effects is the strong dependence of the ion- 
ization (both for projectile and target) cross section, on 
the orientation of the molecular axis with respect to the 
direction of the collision velocity. Earlier the existence 
of such an effect was mentioned in the Refs. @-[||, where 
the ionization of molecular targets in their collisions fast 
highly charged ions is studied. 

Another important effect that appears in the collision 
of fast highly charged projectiles with molecules is so- 
called collision multiplicity effect It implies that 
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after the collision with the first atom of target before the 
relaxation the projectile collides with the second atom 
being in the excited state. Considerable contributions 
by multiplicity effects to the ionization cross section and 
energy losses of fast highly charged ions with diatomic 
molecules [l(| and nanoparticles [TTj] was found in re- 
cently 

In particular, it was shown in the Refs. [T3, [H| that 
the orientation and multiplicity effects will occur when 
the time interval between two subsequent collisions is 
less than the relaxation time. Furthermore, it is clear 
that the multiple collision of the projectile with target 
molecule atoms occurs when the direction of the projec- 
tile motion is close to the orientation of target molecule 
axis. Qualitatively, the role of the orientation effects can 
be understood for diatomic molecule as follows. There- 
fore in the study of fast collisions of highly charged pro- 
jectiles with diatomic molecules one should take into ac- 
count the two-step processes which include collisions of 
the projectile with the first atom of the molecule and 
cannot relax into the ground state before the collision 
with the second atom of the target. We note that in the 
present paper we consider the case when the time be- 
tween two subsequent collisions is much shorter than the 
characteristic period of the projectile electrons. There- 
fore the processes we deal with can be interpreted as the 
variation of the " atomic-double-slit" processes in which 
projectile electrons interacts with both molecular centers 
in a coherent way (see, e.g., [Hf). It is important to note 
that the collision multiplicity effects contribute only to 
the projectile ionization/excitation cross section without 
causing any changes in the transitions of target states. 
Indeed, in subsequent collisions with each atom of the 
target molecule the ionization or excitation occurs in dif- 
ferent atoms. It is clear that the above arguments are 
true for polyatomic molecules, too. 

In this paper we develop a nonperturbative approach 
which is based on the use of sudden perturbation approx- 
imation and utilizes it for the calculations of the electron 
loss cross sections of fast highly charged projectiles in the 
collisions with polyatomic molecules. The method allows 
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to take into account all the transitions both in target 
and projectile electronic states. Moreover, it is possible 
to achieve considerable simplifications of the expressions 
for the electron loss cross sections in the cases of high 
enough projectile charges and multi-electron targets . 

The calculations of the electron loss cross sections and 
their dependence on the orientation of the molecular axis 
are presented for single-, double-, and triple- stripping of 
the fast highly charged iron ions. It is also shown that 
the collision multiplicity effect leads to considerable dif- 
ference between the stripping cross sections for the cases 
of parallel and perpendicular orientation of the molecu- 
lar axis with respect to the direction of projectile mo- 
tion. However, for chaotic orientation of the axis of tar- 
get molecule such effect is negligible. Finally, it should 
be noted that for the calculation of the multiple strip- 
ping cross section of the fast highly charged projectiles 
in their collisions with multi-electron targets one should 
use non-perturbative methods, since the Born approx- 
imation is not applicable for these processes. Our ap- 
proach also can be considered as a nonperturbative one. 
However, our results concerning on the dependence of the 
cross sections on the molecule's axis orientation can be 
obtained within the perturbative approach, too: at least, 
for a single-electron ionization. 

This paper is organized as follows. In the next sec- 
tion we will present detailed description of the approach 
to be utilized, derivation of the nonperturbative expres- 
sions for the stripping cross section and its applications 
for the calculations of the multiple stripping cross sec- 
tions of fast highly charged projectiles in the collisions 
with the nitrogen molecule, N2 . The section III presents 
discussions of the results, while in the section IV we give 
some concluding remarks. 



II. THEORETICAL BACKGROUND 

The system we are going to treat includes a fast highly 
charged projectile colliding with a molecule that consists 
of multi-electron (number of electrons, Na ^> 1) atoms. 
For such collision system the cross section for excitation 
or ionization of the projectile from state |0) to |fc) (pro- 
vided arbitrary transition can occur in the target state) 
with high accuracy (error is proportional to 1/Na) in the 
Glauber approximation can be written as [f| lioj 
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Here v is the projectile velocity and the x— axis is di- 
rected along the vector v. If b is the impact parame- 
ter, and target position is fixed with one of the nuclei 
being at the origin of coordinate system, then the posi- 
tion of the projectiles nucleus can be given by the vector 
R = (X, b) . Atomic units are used here and in the fol- 
lowing. In Eq. ((T|) the potential U describes interaction 



of the projectile electrons with the target which is con- 
sidered as a lengthy object. In other words, the cross 
section in Eq. (fTJ) is expressed in terms of the charge 
density of target electrons [ISIil). Formally, Eq. © 
has the same form as that for the ionization cross section 
in the "frozen" approximation for the target electrons. 
However, as it was shown in the Ref. [4], for the multi- 
electron targets (Na ^> 1) this formula describes (with 
the error order of X/Na) the transition cross section for 
the projectile electrons summed over the complete set of 
final states of the target electrons. For a molecule consist- 
ing of multielectron atoms the electronic density is almost 
the same as the sum of that for isolated atoms. Therefore 
for such case we can consider the target as consisting of 
isolated and non-interacting atoms [l0|, HH). We describe 
the charge density of the each atom in the target within 
the Hartree-Fock-Slater model [l3| in which the spatial 
charge density can be written as 



3 



(2) 

where Z m is the charge of the m-th atomic nucleus, A OT) j 
and a mj i are the tabulated constants that can be found 
in the Ref. [HI . The potential created by molecule at the 
point r can be written in terms of potentials created by 
each atom of the molecule: 



N Z 

m=l m 



(3) 



where N is the number of atoms in target, d m is the 
distance from the m the nucleus in the target to the ob- 
servation point, r and $ m (r) is the screening function for 
the m-th atom. In the Hartree-Fock-Slater model screen- 
ing function can be written as: 
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Furthermore, let us introduce the following notations: 
r p is the coordinate of the projectile electrons with re- 
spect to its nucleus, (p = 1, 2, Np), Np is the total 
number of projectile electrons, R m = (X m ,h m ) are the 
distances between the projectile nucleus and the nucleus 
of the m-th atom of the target. If b m is the impact pa- 
rameter with respect to m-th atom, then d m — |R m + r p |. 

Therefore the interaction potential between the target 
and projectile electrons can be written as 
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Taking into account Eq. ^ for the "eikonal phase" 



3 



involved into Eq. (Q} we have 
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N 2Z Np 3 
= X] — ~ ^2 ^2 A m,iK Q (a m ,i \b m + s p \) , (5) 
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with Kq (z) being the lowest-order McDonald function 
and s p is the projection of the vector r p onto the im- 
pact parameter plane. As it was shown in the Ref. [6j 
Eq. ([T]) with the " eikonal phase" given in the formula ([5]) 
is applicable for relativistic collisions, too. Thus Eq. ^ 
describes the cross section for the transition of the pro- 
jectile electrons from the initial state [0} to a state |fc), 
under the assumption that arbitrary transitions can oc- 
cur in the states of target electrons. The accuracy of 
the formula depends on the number of target electrons, 
Na i.e., the error is of order of the quantity ~ 1/Na for 
N A > 1. 

In the following we will consider the highly charged 
projectiles whose " visible" (effective) charge, Zp is much 
larger than unity (for instance, for the Fe 10+ projectile 
Zp = 10 while nucleus charge is Z = 26). Then the char- 
acteristic size of the projectile is much less than that of 
an atom of the target molecule. Therefore we can assume 
that the mean field created by atom acts uniformly to the 
projectile electrons that corresponds to expansion of the 
"eikonal phase" in Eq. (pjj in terms of a small parameter, 
Sp/b. This allows us to rewrite Eq. (|TJ) for orthogonal |0) 
and \k) in the following form 
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has the meaning of the momentum transfer to each pro- 
jectile electron due to the collisions with m-th atom of 
the target, K\[z) is the first-order McDonald function. 

Furthermore, we note that Eq. ^ can be applied for 
the collisions of highly charged (Zp ;» 1) projectiles with 
the neutral molecules consisting of multielectron atoms. 
This requires fulfilling by Z t the condition Z t 3> 1, with 
t = 1,2,...,N being the number of atoms in target, i.e., 
Z\ is the nucleus charge of the first atom of the target, 
Zm is the nucleus charge of the iV-th atom etc. 

Since the expression for the cross section is derived 
within the Glauber approximation, the energy of projec- 
tile, E should be much larger than that of projectile- 
target interaction, U, i.e., E 3> U and E » 1, where 
k is the projectile momentum and L is the interaction 
radius for the potential U. In the case of neutral molec- 
ular target whose size is much larger than that of pro- 
jectile, as the quantity L one can take characteristic size 



of target. It is clear that this condition is fulfilled for 
fast highly charged projectile. In addition, to make use 
our approach we have to assume that the sudden pertur- 
bation approximation is valid, i.e., for each target atom 
(interaction) collision time, t c ~ L/v between the projec- 
tile and target is much shorter than the period of most 
fastest (inner) electron, r e , to be ionized: 

t c <C r e . 

Fulfilling of this condition means that target elec- 
trons cannot change their position during the collision 
time. In this case target electrons can considered as 
in fixed positions during the collision For rela- 

tivistic projectiles the above condition can be written 
as \J\ — v 2 /c 2 L/v <C T e . The quantity r e can be esti- 
mated for each fixed collision system. For multielectron, 
collision system, with most of the electrons to be ionized 
being in outer shells we have r e ~ 1. 

As it was mentioned above in order to make use Eq. © 
for the calculation of the projectile stripping cross sec- 
tions in the collision with neutral molecular targets the 
following condition should be fulfilled 

Zp = Z - N P > 1, 

where Np is the number of projectile electrons before the 
collision, Z is the charge of the projectile nucleus. This 
means that (unscreened) charge of projectile should be 
large enough to consider the projectile as highly charged 
particle. Important point in the calculations of projec- 
tile's multiple stripping cross section in the collisions with 
polyatomic molecules is the effects of collision multiplic- 
ity and effects of molecular axis orientations. To include 
these effects into consideration we need to assume that 
two or more atoms in the target are on the same line 
which is parallel to the vector v. To demonstrate col- 
lision multiplicity and molecular axis orientation effects 
we consider simplest target, diatomic molecule. We will 
use notation L to denote the orientation vector is of the 
molecular axis and assume that it is directed along the 
line connecting nuclei of two atoms in the molecule. Then 
the cross section presented by formula ^ is the function 
of vector L i.e., a = cr(L) = <j{6,<j>). Furthermore, we 
represent vector L in terms of spherical coordinates with 
the angles <fi, 9 and assume the molecular axis to be di- 
rected along the vector v. 

In the following we will be interested in the stripping 
cross section averaged over the angle (f>: 



<0) = JL / a(9, 



(8) 



Angle will be called the orientation angle of the 
molecular axis. Then the multiplicity effect can be char- 
acterized in terms of 8— dependent relative correction, 8 
defined as 5 — (<?(0) — a±)/a±, where a± is the strip- 
ping cross section at 8 = ir/2. It is easy to see that 
a± described the collision when target moleculae axis is 
perpendicular to the projectile direction. 



4 



III. RESULTS AND DISCUSSION 

Let us now apply Eq. © to the collision system con- 
sisting of fast Fe q+ projectile and _/V 2 molecule. For sim- 
plicity we consider the cases when q — 24, 25, 26. Let us 
start from the most simplest case: stripping of hydrogen- 
like projectile in the collision with diatomic molecule for 
which the cross section can be written as 

<r(6,(f>) = J d 2 b J dk I (k I exp{-z( qi +q 2 )r} | 0) | 2 . 

(9) 

Here r is the coordinate of the projectile electron with 
respect to the projectile nucleus, k is the momentum of 
the projectile electron lost is the momentum transfer 

that can be written as q^ = ^ A J2'i=i a i^i^i( a i^j)T'^ 
where hj is the impact parameter with respect to jih 
nucleus of the target (j — 1,2). 

Using Eqs. ((HJ and ((5} we have calculated the stripping 
cross section a (6), of the hydrogen-like projectile Fe 25+ 
in the collision with nitrogen molecule N2, for different 
orientation angles and collision energies. A quantity we 
are interested to analyze is the correction to the projectile 
stripping cross section due to the collision multiplicity 
which is given as 6(9) = (cr(0) — cr±)/a±. In Fig. 1 the 
results of calculation of such quantity for collision system 
Fe 25+ — N2 are presented. Continuous line in this figure 
represents 6(9) for the collision energy 10 MeV/u, while 
long-dashed and short-dashed lines are the results for the 
energies 100 MeV/u and 1000 MeV/u, respectively. The 
angle 9 is given in radians. 

1 I 1 — 1 — 1 — 1 1 1 1 1 1 1 — 1 — 1 — 1 1 1 1 1 1 1 — 1 — 1 — 1 1 1 1 1 1 — 
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FIG. 1: The dependence of the quantity 6(9) = (a(9) - 
a±)/a± on the orientation o target molecule axis 9 for the col- 
lisions system Fe 25+ —¥ N2 with a± = a(9 = 7r/2). Continu- 
ous line is the result for the collision energy 10 MeV/u; long- 
dashed line is the result for the collision energy 100 MeV/u 
and short-dashed line is 6(9) for the energy 1000 MeV/u. The 
values of 9 are given in radians. 

Now let us consider stripping of helium-like projectile 



Fe 25+ — N2, in the collision with N2 molecule. In this case 
according to Eq. © the cross section for the transition 
of projectile from state |0,0) to \n\,n2) can be written 
as 

a(6, <p) = J d 2 b I (m, n 2 I e-^+^ir^) | Qj ^ p ^ 

(10) 

with ri and r 2 being the coordinates of the projectile 
electrons with respect toprojectile nucleus. 

Following the Refs. [H, [l6|) in our calculations (fi- 
nal and initial) two-electron states of the projectile 
are described as a symmetric products of one-electron 
hydrogen-like wave functions with effective charges equal 
to degree of ionization. Single a 1+ (9), and double a 2+ (9), 
electron loss cross sections of can be calculated using 
Eq. (jlOp . Here the cross section for double stripping 
a 2+ (6), corresponds to transition of both electron into 
the continuum state, while the single projectile ioniza- 
tion cross section a 1+ (9), describes transition of one of 
the electrons into the continuum, by exciting another 
one into a state of the discrete spectrum. Correspond- 
ingly, the quantity 6 can be estimated for both cases. In 
Figs. 2a, 2b, and 2c we plotted 6 as a function of orienta- 
tion 6 for the collision energies 10, 100, and 1000 MeV/u, 
respectively. Continuous line in this figure is the result 
for single electron loss, 6 = (a 1+ (9) — <t_ l + )/<7_ l + , while 
dashed line describes calculation of 6 for double ioniza- 
tion, i.e., 6 = (a 2 +(9) - a 2 +)/a 2 +. 

Finally, we utilized Eq. ((5]) for the calculation of mul- 
tiple stripping of three-electron (Lithium-like) projectile, 
Fe 23+ in the collision with N 2 molecule. In this case the 
loss cross section describing transition from state |0, 0, 0) 
to |m, n2,fi3) can be written as 

= J d 2 b\ (ni, n 2 , n 3 | e - i ^ + ^ ri+l2+r3) | 0, 0, 0) | 2 , 

(11) 

where ri, r 2 , and r3 are coordinates of projectile elec- 
trons. 

As in the case of helium-like projectile, the wave func- 
tions of three electron state are taken as the symmet- 
ric product of one-electron (hydrogen-like) wave func- 
tions with effective charges equal to degree of ioniza- 
tion discussed in the Refs. [H, EH). Single, double, 
and triple electron loss cross sections of Fe 23+ projec- 
tile can be calculated using Eq. (fTTI) . Single electron 
loss implies ionization of one electron, while other two 
electrons can be excited into any state of the discrete 
spectrum. Similarly, in the case if double electron loss 
two electrons are lost, while other electron excited into 
any bound state. As in the cases of single and double 
electron losses, we have calculated the quantity, 6 for 
the cases of single-, double-, and triple-electron losses. 
Figs. 3a, 3b, and 3c present the results of calculation of 
6(9) for the collision energies 10, 100, and 1000 MeV/u, 
respectively. Continuous line in this figure is the result 
for single electron loss, 6 — (a 1+ (9) — <7_ L + )/a-_ L + , long 
dashed line describes calculation of 5 for double ioniza- 
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tion, 5 = (a + (9) — OjJ ) / c± an d short-dashed line is 
the plot of 5 = (a 3+ (9) - a 3 +)/a 3 +. 
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FIG. 2: The dependence of the quantity 5(6) = (a(0) - 
a±)/a± on the orientation o target molecule axis 9 for the 
collision system Fe 24+ — > N2 for the energies 10 (Fig. 2a), 
100 (Fig. 2b), and 1000 MeV/u (Fig. 2c). Continuous line is 
the result for single electron loss, dashed line is the result for 
double electron loss by projectile. 
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FIG. 3: The dependence of the quantity 5(6) = (a(6) - 
a±)/<j± on the orientation o target molecule axis 9 for the 
collision system Fe 2i+ — > N2 for the energies 10 (Fig. 3a), 
100 (Fig. 3b), and 1000 MeV/u (Fig. 3c). Continuous line 
is the result for single electron loss, long-dashed line is the 
result for double electron loss and short-dashed line for triple 
electron loss by projectile. 

As is seen from the Figs. 1-3, for all collision ener- 
gies stripping cross sections considerably depend on the 
orientation angle 9, the difference between "perpendicu- 
lar" and "parallel" orientations is about 50%-100%. It is 
easy to estimate from the geometrical analysis the value 
of the angle 9 at which multiplicity effect becomes essen- 
tial. Since in all cases projectile collides with the neutral 
atoms of the target molecule whose sizes ~ 1, for the 
internuclear distance of target atoms denoted by L, the 
orientation can be estimated as 9 < l/L. Therefore for 
the nitrogen molecule, for instance, we have L = 2.07, 
that gives the estimate 9 < 0.5. This can be seen also 
from the Figs. 1-3. 
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Also, it should be noticed that the behavior of the 
quantity 5 almost the same for all collisions partners. 
To check this we have done (similar to the above) cal- 
culations for one-, two-, and three-electron Ni and Xe 
projectiles colliding with N 2 , O2, and Au 2 targets. The 
difference from the above treated collision systems we ob- 
served for these systems was too small to present them 
in this paper. 

Finally, since in the experiment the cross section for 
chaotic orientation are usually measured, we have done 
calculations of the stripping cross section averaged over 
the angle 8 assuming uniform distribution on the orien- 
tation angle: 

7T 

= J a(L)^= J a(6)^sm6d6. 


As showed the results of such calculations the cor- 
rections due to the collision multiplicity for the case of 
chaotic orientation is too small: order of 0.1 percents. In 
this case the difference between the cross sections a and 
(7 j_ is very small. Thus the collision multiplicity effect 
is considerable only for the case of regular (non-chaotic) 
orientation of the molecular axis, while for chaotic orien- 
tation it becomes negligible. 

IV. CONCLUSIONS 

The electron losses of fast highly charged projectiles 
in the collisions with neutral molecules has been studied. 



Based on the Glauber approximation a nonperturbativc 
approach is developed to estimate single and multiple 
stripping cross sections. Using the method single, double, 
and triple electron loss cross sections of the fast Fe 25+ , 
Fe 24+ , and Fe 23+ ions in the collisions with N 2 molecule 
are calculated. The effect of collision multiplicity, caused 
by the collisions of the projectile with separate atoms of 
the target is analyzed. It is shown that multiplicity ef- 
fect is essential for the case when the orientation of the 
target molecule axis is parallel or perpendicular to the 
projectile velocity direction, while for chaotic orientation 
such effect is negligible. We note that in all the cases 
we calculate the total electron loss cross sections, i.e., we 
perform summation (integration) over the all continuum 
states of the electron after the ionization. Therefore, in 
the present paper we didn't discuss the energies and di- 
rections of the emitted electrons. The above developed 
method is rather simple for application and can be used 
for any (polyatomic) target molecules and for the projec- 
tiles of arbitrary high (including relativistic) velocities. 
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